00 

o 

O 



NLO QCD corrections to pp -^ WW + jet + X 

Stefan Dittmaier", Stefan Kallweit" and Peter Uwer'' * 

° Max-Planck- Institut fiir Physik (Werner-Heisenberg-Institut), D-80805 Miinchen, Germany 

^ Institut fiir Theoretische Teilchenphysik, Universitat Karlsruhe, D-76128 Karlsruhe, Germany 






> 

00 

o 
p 

o 

oo 

O 



X 



We report on a calculation of the next-to-leading order QCD corrections to the production of W-boson pairs 
in association with a hard jet at hadron colliders, which is an important source of background for Higgs and 
new-physics searches at the LHC. If a veto against the emission of a second hard jet is applied, the corrections 
stabilize the leading-order prediction for the cross section considerably. 



1. Introduction 

The search for new-physics particles — 
including the Standard Model Higgs boson — will 
be the primary task in high-energy physics after 
the start of the LHC that is planned for 2008. 
The extremely complicated hadron collider envi- 
ronment does not only require sufficiently precise 
predictions for new-physics signals, but also for 
many complicated background reactions that 
cannot entirely be measured from data. Among 
such background processes, several involve three, 
four, or even more particles in the final state, 
rendering the necessary next-to-leading-order 
(NLO) calculations in QCD very complicated. 
This problem lead to the creation of an "experi- 
menters' wishlist for NLO calculations" |1I2J that 
are still missing for successful LHC analyses. The 
process pp — » W+W~-|-jet-t-X made it to the top 
of this list. 

The process of WW-|-jet production is an im- 
portant source for background to the produc- 
tion of a Higgs boson that subsequently decays 
into a W-boson pair, where additional jet activ- 
ity might arise from the production or a hadron- 
ically decaying W boson ^. WW-fjet produc- 
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tion delivers also potential background to new- 
physics searches, such as supersymmetric parti- 
cles, because of leptons and missing transverse 
momentum from the W decays. Besides the 
process is interesting in its own right, since W- 
pair production processes enable a direct pre- 
cise analysis of the non-abelian gauge-boson self- 
interactions, and a large fraction of W pairs will 
show up with additional jet activity at the LHC. 
Last but not least WW+jet also delivers the 
real-virtual contributions to the next-to-next-to- 
leading-order (NNLO) calculation of W-pair pro- 
duction, for which further building blocks are pre- 
sented in Ref . (3| ■ 

In these proceedings we briefly report on our 
recent calculation [5: of NLO QCD corrections 
to WW-fjet production at the Tevatron and the 
LHC, but here we discuss results for the LHC 
only. Parallel to our work, another NLO study [3] 
of pp -^ W+W"+jet+X at the LHC appeared. 

A tuned comparison of our results with results 
of Campbell et al. [6J and Binoth et al. [71 is 
in progress. Some details on the status of this 
comparison can be found in Ref. [2]. 

2. Details of the NLO calculation 

At leading order (LO), hadronic WW-|-jet pro- 
duction receives contributions from the partonic 
processes qq — > W+W^g, qg -^ W+W~q, and 
qg -^ W"'"W~q, where q stands for up- or down- 
type quarks. Note that the amplitudes for q = 
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u, d are not the same, even for vanishing light 
quark masses. All three channels are related by 
crossing symmetry. The LO diagrams for a spe- 
cific partonic process are shown in Figure [T] 

In order to prove the correctness of our re- 
sults we have evaluated each ingredient twice 
using independent calculations based — as far as 
possible — on different methods, yielding results 
in mutual agreement. 

2.1. Virtual corrections 

The virtual corrections modify the partonic 
processes that are already present at LO. At 
NLO these corrections are induced by self-energy, 
vertex, box (4-point), and pentagon (5-point) cor- 
rections. For illustration the pentagon graphs, 
which are the most complicated diagrams, are 
shown in Figure [2] for a specific partonic chan- 
nel. At one loop WW-|-jet production also serves 
as an off-shell continuation of the loop-induced 
process of Higgs-|-jet production with the Higgs 
boson decaying into a W-boson pair. In this 
subprocess the off-shell Higgs boson is coupled 
via a heavy-quark loop to two giuons. 

Version 1 of the virtual corrections is essen- 
tially obtained as for the related processes of ttH 
[8] and tt-|-jet [9] production. Feynman diagrams 
and amplitudes are generated with FeynAits 1.0 
[To] and further processed with in-house Math- 
ematica routines, which automatically create an 
output in Fortran. The IR (soft and coUinear) 
singularities are treated in dimensional regular- 
ization and analytically separated from the finite 
remainder as described in Refs. |8|llj . The pen- 
tagon tensor integrals are directly reduced to box 
integrals following Ref. [T^. This method does 
not introduce inverse Gram determinants in this 
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Figure 1. LO diagrams for the partonic process 
uu -^ W+W-g. 



Figure 2. Pentagon diagrams for the partonic 
process uii — > W+W~g. 



step, thereby avoiding numerical instabilities in 
regions where these determinants become small. 
Box and lower-point integrals are reduced a la 
Passarino-Veltman [T3] to scalar integrals, which 
are either calculated analytically or using the 
results of Refs. [Mj. Sufficient numerical stabil- 
ity is already achieved in this way, but further 
improvements with the methods of Ref. [TS] are 
in progress. 

Version 2 of the evaluation of loop diagrams 
starts with the generation of diagrams and am- 
plitudes via FejnArte 3.2 [TH], which is inde- 
pendent of version 1.0 [TU]. The amplitudes 
are further manipulated with FormCalc 5.2 [T7] 
and eventually automatically translated into For- 
tran code. The whole reduction of tensor to 
scalar integrals is done with the help of the 
LoopTools library [T7], which also employs the 
method of Ref. [H] for the 5-point tensor inte- 
grals, Passarino-Veltman [T3| reduction for the 
lower-point tensors, and the FF package [18] for 
the evaluation of regular scalar integrals. The di- 
mensionally regularized soft or collinear singular 
3- and 4-point integrals had to be added to this li- 
brary. To this end, the explicit results of Ref. [TT] 
for the vertex and of Ref. [19] for the box inte- 
grals (with appropriate analytical continuations) 
are taken. 

2.2. Real corrections 

The matrix elements for the real correc- 
tions are given by ^ W+W^qqgg and 
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process 



o'ipb] q-shcrpaipb] Acr/stat. error 



pp^WW+ljet + X 46.453(16) 46.4399(94) +0.70 

pp ^ WW + 2jets + X 31.555(17) 31.5747(63) -1.08 



Table 1 

Coniparison of LO cross sections with Sherpa (taken from Ref. [27]). 



-^ W+W^qqq'q' with a large variety of 
flavour insertions for the light quarks q and 
q'. The partonic processes are obtained from 
these matrix elements by all possible crossings 
of quarks and gluons into the initial state. 
The evaluation of the real-emission amplitudes 
is performed in two independent ways. Both 
evaluations employ (independent implemen- 
tations of) the dipole subtraction formalism 
[20| for the extraction of IR singularities and 
for their combination with the virtual corrections. 

Version 1 employs the Weyl-van-der-Waerden 
formalism (as described in Ref. "^) for the cal- 
culation of the helicity amplitudes. The phase- 
space integration is performed by a multi-channel 
Monte Carlo integrator [2? with weight opti- 
mization |23j written in C++, which is con- 
structed similar to RacoonWW [24]. The results 
for cross sections with two resolved hard jets 
have been checked against results obtained with 
Whizard 1.50 [25] and Sherpa 1.0.8 [26]. Details 
on this part of the calculation can be found in 
Ref. [27], the comparison to Sherpa results is 
briefly illustrated in Table \Tr\ In order to im- 
prove the integration, additional channels are 
included for the integration of the difference of 
the real-emission matrix elements and the sub- 
traction terms. 

Version 2 is based on scattering amplitudes 
calculated with Madgraph [28] generated code. 
The code has been modified to allow for a non- 
diagonal quark mixing matrix and the extraction 



^The input parameters of Ref. 1271 are set as below, 
apart from Qs{Mz) = 0.1187 (l-loop evolved to the scale 
/iren = Mfact = A^w)j and a CKM matrix in Wolfen- 
stein parametrization (to 2nd order in A) with A = 0.22. 
The transverse momenta of additional jets are restricted 
by PT.jet > 20GeV, and the jet-jet invariant mass by 
M(jet,jet) > 20GeV. No jet algorithm is applied, since 
genuine LO quantities are considered in Ref. 1271 . 



of the required colour and spin structures. The 
latter enter the evaluation of the dipoles in the 
Catani-Seymour subtraction method. The eval- 
uation of the individual dipoles was performed 
using a C++ library developed during the cal- 
culation of the NLO corrections for tt-|-jet [S]. 
For the phase-space integration a simple mapping 
has been used where the phase space is generated 
from a sequential splitting. 

3. Numerical results 

We consistently use the CTEQ6 [53] set of par- 
ton distribution functions (PDFs), i.e. we take 
CTEQ6L1 PDFs with a l-loop running as in 
LO and CTEQ6M PDFs with a 2-loop running 
as in NLO. We do not include bottom quarks 
in the initial or final states, because the bottom 
PDF is suppressed w.r.t. to the others; outgo- 
ing bb pairs add little to the cross sectior[fl and 
can be experimentally further excluded by anti- 
b-tagging. Quark mixing between the first two 
generations is introduced via a Cabibbo angle 
9c = 0.227. In the strong coupling constant the 
number of active flavours is A^f = 5, and the 
respective QCD parameters are Ag^ = 165 MeV 
and A|^^ = 226 MeV. The top-quark loop in the 
gluon self-energy is subtracted at zero momen- 
tum. The running of as is, thus, generated solely 
by the contributions of the light quark and gluon 
loops. The top-quark mass is mt — 174.3 GeV, 
the masses of all other quarks are neglected. 
The weak boson masses are Afw = 80.425 GeV, 
Mz = 91.1876 GeV, and A/h = 150 GeV. The 
weak mixing angle is set to its on-shell value, i.e. 
fixed by c^ = 1 - s^ = M^/M^, and the elec- 
tromagnetic coupling constant a is derived from 



•^Sizeable contributions result from top-quark resonances 
in the subproccsses pp -^ W~t-|-X, W+t + X,tt -|-X with 



subsequent top-quark decays t — » W+b, t ^ W b, 
are usually treated as separate classes of processes. 
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Figure 3. LO and NLO cross sections for WW+jet production at the LHC: scale dependence with 
renormalization and factorization scales set to n for pTjot.cut = 50 GeV and 100 GeV (taken from Ref. [S]). 



Fermi's constant G^ = L16637 x 10"^ GeV"^ 
according to a = V2GpA/^s^/7r. 

We apply the jet algorithm of Ref. [30] with 
_R = 1 for the definition of the tagged hard jet and 
restrict the transverse momentum of the hardest 
jet by PTjet > PT,jct,cut- In contrast to the real 
corrections the LO prediction and the virtual cor- 
rections are not influenced by the jet algorithm. 
In our default setup, a possible second hard jet 
(originating from the real corrections) does not 
affect the event selection, but alternatively we 
also consider mere WW+jet events with "no 2'^'^ 
separable jet" where only the first hard jet is 
allowed to pass the PTjet cut but not the second. 

Figure [3] shows the scale dependence of the in- 
tegrated LO and NLO cross sections at the LHC 
for pTjct,cut — 50 GeV and 100 GeV. The renor- 
malization and factorization scales are identified 
here {p = /iron = /^tact), and the variation ranges 
from /i — Afw/10 to /^ = 10 Mw- The depen- 
dence is rather large in LO, illustrating the well- 
known fact that the LO predictions can only pro- 
vide a rough estimate. Varying the scales simul- 
taneously by a factor of 4 (10) changes the LO 
cross section by about 35% (70%). 

Only a modest reduction of the scale depen- 
dence to 25% (60%) is observed in the transi- 
tion from LO to NLO if W pairs in association 



with two hard jets are taken into account. This 
large residual scale dependence in NLO, which is 
mainly due to qg-scattering channels, can be sig- 
nificantly suppressed upon applying the veto of 
having "no 2"^^ separable jet" . In this case the 
uncertainty is 10% (15%) if the scale is varied by 
a factor of 4 (10). The relevance of a jet veto in 
order to suppress the scale dependence at NLO 
was also realized [3T] for genuine W-pair produc- 
tion at hadron colliders. 

Further on, we show the integrated LO and 
NLO cross sections as functions of PTjGt,cut in 
Figure m The widths of the bands, which cor- 
respond to scale variations within Mw/2 < /i < 
2Af-vv, reflect the behaviour discussed above for 
fixed value of pT,jet,cut- For the LHC the reduc- 
tion of the scale uncertainty is only mild unless 
WW-|-2jets events are vetoed. 

Finally, Figure [5] shows the pTjet-distribution 
for the differtial LO and NLO cross sections again 
both for pT,jet,cut = 50 GeV and 100 GeV. Note 
that in the two plots the LO and the more in- 
clusively defined NLO distributions are the same 
up to numerical fluctuations, whereas the more 
exclusive predictions differ in the two plots, since 
the veto applied on a second jet — which is not 
present for the two other curves — depends on the 
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Figure 4. LO and NLO cross sections for WW+jet production at the LHC: dependence on pxjot.cut 
(taken from Rcf. ^Sj). 
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Figure 5. px.jct-distribution in LO and NLO for WW+jet production at the LHC: Here again we set 
renormahzation and factorization scales equal to fj, = Mw- The lower plot shows the iiT-factor for both 
definitions of the NLO observables. The curves for the LO and the more inclusive NLO cross section 
("incl") agree in the pTjct-region covered by both plots, whereas the more exclusively defined NLO cross 
sections ("excl") differ due to the definition of the observable. 
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paring the PTjct.cut plot of Figure 3] with a cor- 
responding plot calculated by summing over the 
distributions of Figure O leads to agreement only 
in the observables with no veto on a second hard 
jet applied. 
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